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Abstract 

It is shown that at finite temperature and in the presence of magnetic sources 
magnetic fields are screened. This is proven within the framework of classical 
transport theory both for the Abelian and non-Abelian plasmas. Magnetic 
screening arises in this formalism as a consequence of polarization effects 
occurring in the plasmas, and it is proportional to the inverse of the gauge 
coupling constant. It is then discussed whether this mechanism could be 
relevant in realistic quantum gauge field theories, such as QCD. 
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I. INTRODUCTION 



The studies of gauge field theories at high temperature T have attracted much atten- 
tion in the recent past JI|. The behavior of certain gauge theories, such as QCD, changes 
dramatically in their low or high temperature regimes. Thus, it is generally believed that 
QCD is in a unconfined phase above a critical temperature T c ~ 200 MeV. In principle, due 
to asymptotic freedom, perturbation theory in the high T regime of QCD could be naively 
expected to be valid, but this turns out not to be the case. 

In the high T regime of QCD, thermal excitations produce a plasma of charged par- 
ticles which screens color electric fields. The electric thermal mass can be computed at 
one-loop order in perturbation theory, and it is proportional to gT, where g is the gauge 
coupling constant. A self-consistent inclusion of the color electric mass in the Feynman loop 
computations requires the use of resummed perturbation theory 0. 

Color magnetic fields are not screened at the same order of perturbation theory. Due to 
this fact several infrared (IR) divergences are encountered in the computations of different 
physical quantities. Those divergences even appear in the absence of quark matter in the 
computation of the perturbative partition function Q. 

There has been a lot of discussions in the literature about the possible mechanism that 
could cure the IR problems in the magnetic sector of a pure Yang-Mills theory. It is generally 
believed @ that non-perturbative effects generate a color magnetic mass of order g 2 T, which 
would arise only in the non-Abelian theory but not in the Abelian case. The arguments to 
reach to that conclusion are the following. The infrared limit of QCD at high T is governed 
by the spatial vector gauge field A^. Using imaginary time formalism, the high T limit of 
QCD is equivalent to an Euclidean three dimensional Yang-Mills theory, with an effective 
coupling constant g' 2 = g 2 T. Those theories are believed to generate a dynamical mass gap 
proportional to the (dimensional) coupling constant, which is not computable in perturbation 
theory. 

The above arguments are rather qualitative, but give no clue about what kind of mech- 
anism could be responsible for the generation of the magnetic mass gap. They also do not 
explain how those effects could eliminate the IR problems found in specific computations. 

Some attempts have been made in the literature to give an effective action that describes 
a magnetic mass term and such that could be used in specific Feynman loop diagrams 
computations However, no conclusive evidence about the proposed Ansatze has been 
reached. 

On the other hand, some discussions have been raised in the literature about whether the 
high T limit of QCD could be correctly described within the classic rather than quantum 
field theory. In Ref. || it has been realized that both the color electric and magnetic 
screenings are detected in lattice classical gauge field computations. This fact could be a 
clear indication that the magnetic screening could be understood in terms of classical but 
non-perturbative physics. 

The purpose of this article is to investigate under which circumstances static magnetic 
screening can emerge in a hot non-Abelian theory using exclusively classical (or semiclassical) 
physics. I do not attempt here to solve the infrared problem of a Yang-Mills theory at finite 
temperature, but I just try to find a classical mechanism that could generate a thermal color 
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magnetic mass. 

A study of the classical gauge field equations reveals that the absence of color magnetic 
screening in a hot plasma is entirely due to the absence of static magnetic sources. Therefore, 
in order to derive a thermal magnetic screening, I postulate the existence of magnetic charges. 
It will be shown in this article that in the presence of magnetic charges magnetic fields are 
naturally screened at finite T. The mechanism responsible for this screening is actually the 
same one as that generating the screening of electric fields in the QED and QCD plasmas. In 
hot plasmas of electric charges, polarization phenomena screens electric fields. The same kind 
of polarization effects would generate the magnetic screening if magnetic charges existed. 

Classical transport theory will be used to prove that the existence of (non-) Abelian 
magnetic charges at finite T implies the screening of (non-) Abelian magnetic fields. This 
formalism has already been used to derive the screening of (non-) Abelian electric fields, and 
those effects are reproduced exactly in the corresponding quantum field theory at high T. 
The generalized set of (non-) Abelian Vlasov equations in the presence of magnetic sources 
will be written, and from them magnetic screening will be derived. As it will be shown, the 
magnetic screening turns out to be proportional to the inverse gauge coupling constant. 

This paper is structured as follows. In Sec. pi], the Abelian plasmas are first studied. 
In Subsec. II A it is recalled how Debye screening is obtained from the Vlasov equations. 
In Subsec. |II B| the proposed "magnetic" or dual Vlasov equations are written, and from 
them the static screening of magnetic fields is derived. It is stressed there that the duality 
symmetry of electromagnetism allows to derive the magnetic screening from the electric 
one. In Sec.|T| the same study is reproduced for non- Abelian plasmas, when duality is not 
a symmetry of the theory. Subsec. [IV A| is devoted to review some static magnetic monopole 
fields solutions in the vacuum. In Subsec. P"V B| solutions to the gauge field equations in 
hot plasmas, which reproduce the screened magnetic fields, are found. Sec. |V| ends with a 
discussion of the results. Let us finally mention that throughout this paper a system of units 
where h = c = ks = 1 will be used. 



II. THE ULTRARELATIVISTIC ABELIAN PLASMAS 



A. Static Electric Screening in the Plasma of Electric Charges 

In this subsection the derivation of the electric Debye screening effects from classical 
kinetic equations is reviewed ||, |7j. This will teach us how to derive the corresponding 
magnetic screening in the presence of magnetic charges. 

Let us first recall the dynamical evolution of a charged point particle. A particle carrying 
an electric charge e, with mass m, and transversing a worldline y^r), where r is the proper 
time, obeys the equations (neglecting the effects of spin) 

jf(r) , (2.1a) 
eF^{y{r))p u {r) , (2.1b) 
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where the electromagnetic field is evaluated on the particle worldline. In our conventions 
F° l = —E\ and F lj = —e tjk B k , where E l and B l are the electric and magnetic fields, 
respectively. 

In a self-consistent picture the electromagnetic fields obey the Maxwell's equations which 
have as sources the electric currents obtained from each charged particle of the system. Thus 

d u F^(x) = r(x) = E E f(*) . (2-2) 

species helicities 

where (helicity and species indexes are implicit) 

f{x) = ej dr^l 5^(x-y(r)) . (2.3) 

The above current is conserved, d^{x) = 0, as may be checked by using the equations of 
motion ( |2.1|) . This is required as a compatibility condition, as can be easily recognized by 
applying a partial derivative <9 M to (|2.2j ), since F^ u is antisymmetric in their indices. 
The electromagnetic field tensor obeys the Bianchi identity 

d v *F^{x) = , (2.4) 

where the dual field is *F V ^ = \e uwa F pa , and e uppa is the Levi-Civita antisymmetric tensor 
in four dimensions, with e 0123 = 1. In our conventions *F 0t = —B % , and *i™ = e^ k E k . The 
Bianchi identity ( |2.4| ) guarantees the existence of the vector gauge field A^(x), such that 
Ffa, = d^A u — d u A^, which is not unique. The ambiguity in defining the vector gauge field 
is that corresponding to gauge transformations. 

The statistical description of the plasma of charged particles is given by the distribution 
function of its components in their phase-space. In the collisionless case the one-particle 
distribution function f(x,p) of finding a particle in the state (x,p) evolves in time via a 
transport equation 

= . (2.5) 
Using the equations of motion (|2.1|) , it becomes the Boltzmann equation 



if 



d „ , v d 
-ei>(s)- 



f(x,p) = 0. (2.6) 



dx^ dp,, 

In a self-consistent picture, the mean electromagnetic fields obey the Maxwell's equations 

d v F^(x) = J»{x) = J2 E 3*i x )i (2-7) 

species helicities 

d v *F"»{x) = , (2.8) 

where now each particle species electric current is obtained from the corresponding distri- 
bution function as 

r(x) = eJdPp»f(x,p). (2.9) 



3 



The momentum measure in ( |2.9| ) is defined as 



dP=-0^2 6( Po )5(p 2 -m*), (2.10) 

so that it guarantees positivity of the energy and on-shell evolution. 
The equations ( |2.6D -( t2~9"|) are known as Vlasov equations ||. 

The effects of static screening in the plasma of electrons and ions can be deduced from 
the Vlasov equations as follows H, 0. Let us consider a neutral plasma, that is, composed 
by the same number of positive and negative charges. The plasma, initially at equilibrium, 
is disturbed by a weak electromagnetic field. We look for a distribution function of the form 

f(x,p) = f^(p ) + ef^(x,p) + ... , (2.11) 

where /'°'(po) is, up to a normalization constant, the Fermi-Dirac equilibrium function 

MPo) = • (2.12) 

Neglecting second order terms, obeys the following equation 

p^fU(x,p)= P »F fl0 (x)J-f(°\p ) . (2.13) 

Notice that only the electric field enters into the r.h.s. of ( |2.13j ), but not the magnetic field. 
This is actually the reason why there is only static electric screening but not static magnetic 
screening in this approach. 

A total electric current density J^(x,p) is defined such that the total current J M (x) is 
found just by integrating over the momenta p, using the momenta measure ( |2.10| ). The 
induced electric current density J fl (x,p) obeys the equation 

p-dJ"(x,p) = eV/F rt (^/(»)( ft ) . (2.14) 

dp 

In the ultrarelativistic limit, that is, taking the fermion mass m = in ( |2.10| ), the induced 
current in momentum space reads 

nk) = -im% f^—f—v-m, (2-15) 
J Arc v ■ k + ze 

where m? D = e 2 T 2 /3 is the Debye mass squared. Retarded boundary conditions have been 
imposed in (|2.15| ), with the prescription ie. The four vector t> M = (1, v) is the four velocity 
of the particles of the plasma, and in the ultrarelativistic situation considered here, it is 
light-like. The angular integral in (|2.15 ) is defined over all possible directions of the three 
dimensional unit vector v. 

In the static situation J 1 = 0, while the induced electric density is 

T n/ \ of d 3 k k • E(k) ,- w , 

J°(x) = -tm% J — — jAie** . (2.16) 
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The Maxwell's equations which have as source the current (|2.15 ) are also known as Kubo 
equations. In the static limit they read 



V • E (x) = J°(x) , V x B (x) = , (2.17a) 
V-B(x) = 0, VxE(x)=0. (2.17b) 

These equations describe the static screening of electric fields inside the plasma. In the 
static situation the magnetic fields obey the same equations as in the vacuum. In the non- 
static situation, magnetic fields also suffer a dynamical screening in the plasma, but we will 
not be concerned in this article about dynamical effects. 

Let us finally recall that the static screening effects described in this subsection have 
been reproduced in the context of perturbative QED in the high temperature limit. 



B. Static Magnetic Screening in the Plasma of Magnetic Charges 

In this subsection the effects of static screening in a plasma of magnetic charges are 
derived. 

Let us first recall the classical equations of motion of a particle carrying a magnetic 
charge e, with mass m, and transversing a worldline y^{r) || 

p"(r) , (2.18a) 

e*F^(y(T))p u (r) . (2.18b) 

In a self-consistent picture, these are augmented with the field equations 

d v F Vft (x)=0, (2.19) 

d u *F^{x) = J"{x) = E ~^( x ) ■ ( 2 - 2 °) 

species helicities 

The magnetic current is computed for each particle species as 

~f{x) = ejdr 6®(x - y(r)) , (2.21) 

and it is conserved d lx j tl {x) = 0. 

Comparing Eqs. (|2.1J )-( j2.4|) and ( [2.18j )- (p.21| ), we see that they are symmetric under the 
interchange of electric and magnetic fields (E, B) — > (B, — E), and electric charges by mag- 
netic ones. This is the so called duality symmetry of electromagnetism. 

In the presence of a magnetic charge, and due to the absence of the Bianchi identity 
( p.4| ), it is not ensured that the electromagnetic field can be derived globally from a vector 
gauge field as F^ u = d fl A u — d v A fM . However, one can still define a vector gauge field which 
obeys that condition outside the worldline of the magnetic charge. Therefore, it is possible 
to define a vector gauge field locally. For the case of the point magnetic charge, it is 



m 



m ■ 
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dr 
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enough to define two different vector gauge fields in different space-time regions, the two 
solutions being related in their common domain of definition by a gauge transformation^. 
However, it would be possible to derive the dual electromagnetic field from a dual vector 
gauge field A v , such that *F flu = d^A v - d v A^ since dT^ = 0, and = -^ ual3 *F a/3 . 

The above equations show the exact duality of electromagnetism when electric and mag- 
netic degrees of freedom are interchanged. Due to this fact, one naturally expects that the 
propagation properties of electric fields in the plasma of electric charges are the same as 
those of magnetic fields in the plasma of magnetic charges. In particular, this is true for 
static screening. Although from the duality symmetry arguments this is obvious, let us show 
how one would derive the static screening of magnetic fields. 

Let f(x,p) be the probability distribution function of finding a particle with mass m 
carrying a magnetic charge e in the state (x,p). The corresponding set of dual Vlasov 
equations can be derived from the equations of motion ( p.lSQ , and it reads 



pi-'- 



ex „ - -^-^ />.rt = 0, (2.22) 
d u F Ufl {x) = , (2.23) 
d v *F^{x) = J»{x)= J2 E (2-24) 

species helicities 

Now each particle species magnetic current is obtained from the corresponding distribution 
function as 

~f(x) = ej dPp^f(x 1 p) . (2.25) 

Let us consider a neutral plasma of magnetic charges which is initially at equilibrium. If 
the system is perturbed by a weak field, then / can be found in the form 

f(x,p) = / (0) (po) +ep\x,p) + ... , (2.26) 

where is, up to a normalization constant, the Fermi-Dirac equilibrium distribution 
function np(po). The equation obeyed by /W reads 

p^f il \x,p)=p»*F, (x)J-fM(p ) . (2.27) 



Notice that only the magnetic field enters into the r.h.s. of ( 2.27Q , since *F i0 = —B 



In the ultrarelativistic limit, i.e. taking m = 0, the conserved magnetic current in Fourier 
space is given by 



1 Alternatively, one could work with an unique vector gauge field with a Dirac string attached 
to the monopole. We prefer the Wu-Yang H construction which eliminates references to the Dirac 



string (see Subsec. IV A ). 
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J> W = ^./!£_£_ V .B(*), (2.28) 

where m 2 D = e 2 T 2 /3 is the magnetic Debye mass squared. The notation that has been 
used above is the same as in Eq. ( [2.15| ), that is, retarded boundary conditions have been 
implemented, and is the four light-like velocity of the particle. 
In the static situation J 1 = 0, and the dual Kubo equations read 

V • E (x) = , V x B (x) = , (2.29a) 

V • B (x) = J°(x) , VxE(x) = 0, (2.29b) 

where 

J ° (x) = / 'wf^ 1 ^ ■ (2 - 30) 

These equations describe the static screening of magnetic fields. Notice that they could 
have been obtained from Eqs. (|2.17|) just by using duality symmetry arguments. 

Since the plasma is composed of quantum particles, the magnetic charge e is subject to 
the Dirac quantization condition 

ee = 2irn, neZ . (2.31) 

Therefore, the magnetic Debye mass can be written in terms of the electric charge as 

T 

m D = —^2im . (2.32) 
V3 e 

Let us notice that in the perturbative regime e<l, the expansion performed in Eq.( [2.26 ) 
would have failed. 



III. THE ULTRARELATIVISTIC NON-ABELIAN PLASMA 



A. Static Screening of Color Electric Fields in the Plasma of Non-Abelian Charges 

Consider a particle bearing a non-Abelian SU(N) color charge Q a , a = 1,...,N 2 — 1, 
traversing a worldline j/ m (t). The Wong equations [[II]] describe the dynamical evolution of 
the variables x M , p M and Q a (we also neglect here the effect of spin): 

p"(r) , (3.1a) 

gQ a (T)Fr(y(T))p u (r) , (3.1b) 

-SrV(r)4(2/(r))Q c (r) , (3.1c) 



rn 



m ■ 



m 



dy»{T) 
dr 

dr 
d Q a (r) 
dr 
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where f abc are the structure constants of the group, denotes the field strength, which is 
defined as 

F a , v = d,A a u - d v Al + gr hc A\ Al . (3.2) 
The color fields obey the Yang-Mills equation 

[D v F^f{x) = J» a {x) = E E f a ^) > (3-3) 

species helicities 

where the color current associated to each colored particle is 

j;(x)=g J dr Ptl (T)Q a (T)5^(x-y(r)) . (3.4) 

As a consequence of the definition (|3.2| ), the color fields obey the non-Abelian Bianchi 
identity 

[D u *F^] a (x) = , (3.5) 

where *F vpa = \e vpp(J F a pa . 

The main difference between the equations of electromagnetism ( j2.1| ) and the Wong 
equations ( |3.1| ), apart from their intrinsic non-Abelian structure, comes from the fact that 
color charges precess in color space, and therefore they are dynamical variables. Equation 
(|3. lc|) guarantees that the color current associated to each colored particle is covariantly 



conserved 

(D,f) a (x) = d„f a {x) + gf abc Al(x)f c (x) = , (3.6) 
therefore preserving the consistency of the theory. Notice that ( |3.6|) is required as a compat- 



ibility condition, since in applying a covariant derivative to (3^3) the l.h.s of the equation 
should vanish. 

Let us formulate the statistical description of a plasma of colored particles in their phase- 
space. The usual (x,p) phase-space is enlarged to (x,p, Q) by including color degrees of 
freedom for colored particles. Physical constraints are enforced by inserting delta-functions 
in the phase-space volume element dxdPdQ. The momentum measure, chosen as in ( p,10| ), 



guarantees positivity of the energy and on-shell evolution. The color charge measure enforces 
the conservation of the group invariants, e.g., for SU(3), 

dQ = d 8 Q 5(Q a Q a - q 2 ) 5(d abc Q a Q b Q c - q 3 ) , (3.7) 

where the constants q 2 and fix the values of the Casimirs and d abc are the totally symmetric 
group constants. The color charges which now span the phase-space are dependent variables. 



These can be formally related to a set of independent phase-space Darboux variables []TT 
For the sake of simplicity, the standard color charges will be used in the remaining part of 
this Section. 

The one-particle distribution function f(x,p,Q) denotes the probability for finding the 
particle in the state (x,p,Q). In the collisionless case, it evolves in time via a transport 



S 



equation j- = 0. Using the equations of motion ( |3.1| ), it becomes the Boltzmann equation 




r\ r\ r\ 

- g Q a F"(x)—-gf abc A b (x)Q c ' 



dx^ 



dp p 



f(x, P ,Q) = 0. 



(3.8) 



A complete, self-consistent set of non-Abelian Vlasov equations for the distribution func- 
tion and the mean color field is obtained by augmenting the Boltzmann equation with the 
Yang-Mills equations: 



[D v F vfM ] a {x) = J^ a (x) 



E E r a (x) 

species helicities 



(3.9) 



where the color current j^ a {x) for each particle species is computed from the corresponding 
distribution function as 



f a {x)=g JdPdQ p»Q a f(x,p,Q) . 



(3.10) 



Notice that if the particle's trajectory in phase-space would be known exactly, then Eq. 
( P . 1Q ) could be expressed as in Eq. (|3.4| ). Furthermore, the color current ( 3.10|) is covariantly 
conserved, as can be shown by using the Boltzmann equation |Tj|] . 

The Wong equations (|3.1|) are invariant under the finite gauge transformations (in matrix 
notation) 



A, 



uA lt ir 1 --u^-u- 1 

9 ox. 



(3.11) 



Q = UQU' 1 ; 

where U = U(x) is a group element. 

It can be shown [11[ that the Boltzmann equation ( |3.8| ) is invariant under the above 
gauge transformation if the distribution function behaves as a scalar 

f(x,p,Q) = f(x,p,Q) . (3.12) 

To check this statement it is important to note that under a gauge transformation the 



derivatives appearing in the Boltzmann equation (|3.8| ) transform as [11 
2 Tr ([ { ±-U)U-> 



d 



d 



d 



d 



dx^ 1 dx^ 



dpv dpv 



dQ dQ 



(3.13) 



that is, they are not gauge invariant by themselves. Only the specific combination of the 
spacial and color derivatives that appears in (|3.8| ) is gauge invariant. 

The color current ( |3.10| ) transforms under ( |3.11| ) as a gauge covariant vector: 



f{x) 



g J dPdQp^UQU- 1 f(x,p,Q) = U f(x) U' 1 . 



(3.14) 



This is due to the gauge invariance of the phase- space measure and to the transformation 
properties of /. 
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Let us consider now a gluon plasma, that is, a plasma of particles carrying a non-Abelian 
charge in the adjoint representation, and which is initially at equilibrium. The system is 
disturbed by a weak color electromagnetic field, and one looks for the response of the plasma. 
The distribution function can be expanded in powers of g as: 

f = f {0) + gf {1) + ... , (3.15) 

where /(°) is, up to a normalization constant, the Bose-Einstein equilibrium distribution 
function 

= e Po/T _ I ■ ( 3 - 16 ) 

The Boltzmann equation ( |3.8|) for reduces to [11 



P" ~ 3 A >)Q^) f { ^P,Q)=P"Q a F^(x)^-f^{po) ■ (3.17) 

Notice that a complete linearization of the equation in would break the gauge invariance 
of the transport equation, which is preserved in this approximation. But notice as well that 
this approximation tells us that also carries a (^-dependence. 

One can get the equation that the color current density J£(x,p) obeys by multiplying 
( P-17j ) by pt 1 and Q a and then integrating over the color charges. For gluons in the adjoint 
representation 



dQQaQb = N5 ab , (3.18) 
one finally gets, after summing over helicities, 

[p-D J»(x,p)] a = 2g 2 Np^F? (x)-^f(°\p ) . (3.19) 

apo 

Notice that only the color electric field enters in the r.h.s. of the above equation. Thus, only 
the color electric field is screened in the static situation. 



The induced color current can be expressed in terms of the parallel transporter $ as [13 

//7Q - rca 
— -v" / du<5> ab (x,x-vu)v-E b (x-vu) (3.20) 

where m 2 D = g 2 T 2 N/3 is the Debye mass squared, and $ obeys the equation 
d 

—$ ab (x,x - vu) = g§ ac (x,x- vu) fcbdV^A^x - vu) (3.21) 

with the initial boundary condition § a b(x,x) = 5 a b. The four vector t> M is the velocity 
vector of the particles of the plasma, and in the ultrarelativistic limit is light-like. Retarded 
boundary conditions have also been implemented in Eq. (|3.20| ). 

Alternatively, in momentum space J^ik) may be expressed as an infinite power series in 
the vector gauge field A a Jk) fL4f . 
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In the static situation the color current simplifies drastically and may be expressed as 



J» = -m^ Ag(x) 0. 

The non-Abelian Kubo equations were first derived in Ref. |I4 
read 



(D ■ E)„ = J?(x) 



(D B) 







(BxB) a + gf abc A b E c = , 
(DxE) a -gf abc A b Q B c = . 



In the static limit they 

(3.22) 
(3.23) 



The static non-Abelian Kubo equations have been studied in Ref. [[L5| and they describe 
the static screening of color electric fields. Furthermore, the color electric screening effects 
described by these equations are reproduced in the context of resummed perturbative QCD 
in the high temperature limit M. 



B. Static Screening of Color Magnetic Fields in the Plasma of Non-Abelian 

Magnetic Charges 

In this subsection the screening effects in a plasma of non-Abelian magnetic charges 
are derived. Those particles are the natural non-Abelian analogues of the magnetic Dirac 
monopoles that were studied in Subsec. [11B. 



As opposed to what happened in the Abelian theory, one cannot appeal to duality 
symmetry arguments to describe the dynamics of non-Abelian magnetic charges. It is a 
well-known fact that a pure Yang-Mills theory without matter is not symmetric under the 
interchange of color electric and magnetic degrees of freedom |Tol . This may be understood 



as follows. In the absence of matter, the Yang-Mills equation and the non-Abelian Bianchi 
identity both involve a covariant derivative, but not a "dual" covariant derivative. Thus, 
in terms of the dual field *F£ V the Yang-Mills equation cannot be interpreted as the non- 
Abelian Bianchi identity for *F£ V . In general, it is not possible to find a dual vector gauge 
field which is related to *F^ U as A a ^ is related to Fjf u in (|3.2|). It has been realized in the 



literature that the dual of a Yang-Mills theory is a Freedman-Towsend like theory ]T7 



The fundamental object of that theory is an antisymmetric two index tensor, Fj£ v . One can 
define in that theory a vector gauge field potential which acts as a parallel transport 
for the phases of the charged particles, but which is not related to F£ v as in a Yang-Mills 
theory. 

The equations of motion of a particle of mass m, carrying a non-Abelian magnetic charge 
of SU(N) Q a , with a = 1, N 2 — 1, and transversing a worldline y M (r) have been derived 
in Ref. |fL9| . The derivation was made there using a variational principle in loop space. The 
equations written in terms of the color gauge fields are 

m^^=^(r), (3.24a) 



dr 

m dP^ = ~Qa {T y K , {y{T))pv{r) ^ (324b) 
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m 



d Q a {r) 
dr 



-g 



f abc p"(T)A b (y(r))Q c ( 



T 



These are augmented with the field equations 

[D u F v %(x) = , 



[D v *F»e\ a {x) = J*\ 



X) 



E E # 

species helicities 



X) 



where 



~fJx) = g / drp,{r) Q a (r) 8^\x - y(r)) . 



(3.24c) 

(3.25) 
(3.26) 

(3.27) 



The above current is covariantly conserved [D^ J^] a (x) = 0, which is required by consis- 
tency, as may be realized in applying a covariant derivative to ( |3.26| ). 

Notice that a new coupling constant g appears in the equations as coupling between the 
color electromagnetic fields and the non-Abelian magnetic charges. In principle, this is a 
new variable in the system. However, there is a Dirac quantization condition relating g and 
g. We will come back to this point later on. 



The above equations are the dual of the Wong equations written in Subsec. [Ill A| . Notice 
that they are not symmetric under the interchange of color magnetic and electric degrees 
of freedom, since the same covariant derivative appears in both of them. Some caution 
should be taken at this stage. Exactly as it happened in the Abelian situation, in the 
absence of the non-Abelian Bianchi identity one cannot define globally a vector gauge field 

obeying the field equations. It could be defined locally, outside the monopole worldline. 
It is possible to define different in different space-time regions, the different solutions or 
patches being related in their common domain of definition by a gauge transformation. 

Let us stress that the non-Abelian electric and magnetic charges Q a { T ) an d Q a { T ) live in 
the same group manifold, and obey the same kind of dynamical evolution. Both transform 
under gauge transformations in the same way. 

With the above classical equations of motion it is possible to derive the dual non-Abelian 
Vlasov equations. They read 



P 



° gQ a *FUx)^-gf abc A b (x)& ° 



dxv- 

[D u F v »] a {x) = 



dp u 



dQ a 



f(x,p,Q)=0 



[D v *F^]\x) = r a {x)= £ £ ~f\x) 

species helicities 



where 



f a {x) 



g dPdQp»Q a f(x,p,Q) 



(3.28a) 

(3.28b) 
(3.28c) 



(3.29) 



The color magnetic charge measure dQ is defined as in Eq. (p.7|). 
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The gauge transformation properties of the above equations will not be discussed here, 
as they turn out to be exactly the same as their dual partners. 

Let us derive now the response of the non-Abelian magnetic plasma to a weak field 
disturbance. Let us suppose that the neutral plasma of non-Abelian magnetic charges in 
the adjoint representation is initially at equilibrium. Then one can look for a solution of the 
form 

f = f (0) + ~gf {1) + ... , (3.30) 

where is, up to a normalization constant, the Bose-Einstein equilibrium distribution 
function. Then f^ 1 ' obeys the equation 

P" U^-9f abc Al(x)Q c -^j P(x,p,Q) = P »Q a *F« (x)^-f(°\p ) . (3.31) 

A weak coupling expansion in g, the natural coupling constant of the problem, has been 
performed above. To preserve the gauge symmetry of the above Boltzmann equation, the 
term linear in A a in the l.h.s. of (|3.31|) has to be kept. This tells us that has a dependence 



on the coupling constant g. The two constants g and g are related by a Dirac quantization 
condition that will be discussed later on. 

From Eq. (|3.31| ) one can derive the static screening of color magnetic fields. In order to 
do that, one should follow the same procedure as in the previous subsections. Notice that 
now only the color magnetic field enters in the r.h.s of ( p. 31 ). 



The equation obeyed by the current density is obtained after integrating over the mag- 
netic charges Q a , and summing over helicities, and it reads 

[p-DJ»(x,p)} a = 2~g 2 Np^*F* (x)-^-f(°\po) . (3.32) 

dpo 

The solution of the above equation may be written as 

m( x ) = m 2 D / —^v" / du$ ab (x,x -vu)v -B\x -vu) , (3.33) 
J Atx Jo 



where $ is the parallel transporter, which obeys also Eq. ( |3.21|) . Here fh? D = g 2 T 2 N/3 is 
the magnetic Debye mass squarred. 

The dual Non-Abelian Kubo equation read in the static limit, therefore 

(D • E) a = , (D x B) a + gf abc A b E c = , (3.34) 

(D • B) a = J a °(x) , (D x E) a - gf abc A\ B c = , (3.35) 

and they describe the static screening of color magnetic fields. 

There is a Dirac quantization condition relating the two coupling constant g and g, which 
depends on the Lie group under consideration. For the case considered here, where all the 
matter is in the adjoint representation, then the gauge field theory which is based upon the 
Lie algebra su(N) has as a global Lie group SU(N)/Z N and not SU(N). Here Z N is the 
iV-element finite group consisting in the iV-th roots of unity, that is, the integral powers of 
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exp (2-ki/N). This is due to the fact that two SU(N) matrices that differ only by a factor 
belonging to Zn will be represented by the same matrix in the adjoint representation. For 
a SU(N)/Zn theory the Dirac quantization condition reads | 22fl , [|i~9[|, |f2~Tf 



99 



2iin 



neZ . 



Therefore, the dual Debye mass is written in terms of g as 

T 



2tt n 



^3.36) 



(3.37) 



Therefore, for small values of N the static magnetic screening effects could not have been 
reproduced by using perturbation theory in g, as then the expansion of Eq. (|3.30|) would 
not hold. However, notice that for N — > oo, both g and g can both be small, and therefore 
perturbative expansions in g and g can both be valid. 



IV. STATIC MAGNETIC MONOPOLE FIELDS 

In this section exact solutions to the dual Kubo equations are found. Those solutions 
describe screened magnetic fields in the plasma of magnetic charges. Some known results on 
how to construct magnetic monopole fields in the vacuum for the Abelian and non-Abelian 
theories are first reviewed. 



A. Monopole Fields in the Vacuum 

The magnetic field created by a magnetic charge e which is at rest at the origin of 
coordinates is given by 

An r 3 



As already explained in Subsec. I fBl it is not possible to construct a global vector gauge 



field A^(x) which generates the above magnetic field. However, it is possible to find a local 
vector gauge field which is defined everywhere except on a "Dirac string". That was the 
original construction of magnetic monopole fields due to Dirac. Wu and Yang || showed 
that it is also possible to construct vector gauge fields without references to Dirac strings as 
follows. Let us denote by R the space region surrounding the magnetic monopole. Dividing 
R into two regions, R + and -R_ defined as (in spherical coordinates (r, #,</>)), 

R + : < 6 < vr/2 + 5, r > 0, < < 2vr , (4.2) 
R_ : tt/2 - 5 < 9 < n, r > 0, < <p< 2tt , (4.3) 

with an overlap region 7t/2 — 5<9<tt/2 + 5, where < 5 < tt/2. In each region one may 
take [| 
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e (1 - cost)) „ ( a a \ 

A+ = — e , 4.4a) 

Airr sin 

. e (1 + cos#) A N 

A- = - — - e . 4.4b 

47rr sm # 



These two vector gauge potentials reproduce the magnetic field (|4.1| ) in their respective 
domain of definition. Furthermore, in the overlap region the two vector gauge fields are 
related by a gauge transformation 

In a quantum description, this gauge transformation implies a change in the wavefunction 
■0 of a particle , ip — > exp (— zee0/27r)-?/>. After requiring the wavefunction be single-valued, 
one arrives at the Dirac quantization condition ee = 27m. 

One can easily generalize the previous results to the non-Abelian case. Let us consider 
first for definiteness the SU{2) Lie group and take as its infinitesimal generators t a = — |cx a , 
where a a are the Pauli matrices. A non-Abelian magnetic monopole field may be constructed 
from the Abelian one by taking the same value (|4.1|) in a particular direction in color space. 
For example, 

B 3 = -f ^ > B x = B 2 = . (4.6) 
The vector gauge fields which reproduce the above color magnetic field are obtained from 



taking in the a = 3 direction in color space the Abelian vector gauge fields ( |4.4| ), while 
= A% = 0. 

In general, for a simple Lie group, the corresponding magnetic monopole fields can be 
constructed from their Abelian counterparts, just by multiplying the Abelian fields by a 
constant matrix Q living in a specific representation of the Lie group. Those are solutions to 
the Yang- Mills equations in the presence of a static point magnetic source J M = 5^q g Q5^ (r). 

Due to the Abelian character of these non-Abelian monopoles, the Dirac quantization 
condition can be easily worked out, and it reads 

exp ( 2 .f?) = 1 . (4.7) 

The above quantization condition is sensitive to the global structure of the Lie group under 
consideration, and it could be different for different Lie groups sharing the same Lie algebra. 
For example, for SU{2) and 5*0(3), Eq. ( ^4.7|) has different implications |J. 

There is a dynamical [21], as well as topological classification [^2| of these non-Abelian 



monopoles. The topological classification of the monopoles associates each class of monopole 
to each different element of the first homotopy group tt\ (G) of the Lie group G under con- 
sideration. A stability analysis of the non-Abelian magnetic monopole fields was performed 



in [20], finding that there is one stable non-Abelian monopole field for each topological class. 
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B. Static Magnetic Monopole Fields in the Hot Plasmas 



In this subsection some exact solutions to the static dual Abelian and non-Abelian Kubo 
field equations are presented. In order to find those solutions, rather than solving the 
corresponding differential equations, we will take profit of the following facts. First, in the 
Abelian situation, the solutions to the "magnetic" Kubo equations can be found from those 
of the "electric" ones, just by making use of the duality symmetry of electromagnetism. 
Second, solutions to the non-Abelian Kubo equations can always be constructed from the 
Abelian ones, as we have already explained. Using these facts, one can easily construct 
solutions to the dual non-Abelian Kubo equations which have an Abelian character. 

Let us first recall how to find exact solutions to the static Kubo equations ( |2.17 ). Those 



can be easily solved in terms of the electric potential A , with E = — VA . Then, the 
equation obeyed by A is 



V 2 - m 2 D ) A (x) = , (4.8) 



With spherical symmetry, and in spherical coordinates (r, #,0), Eq. ( |4.8j ) can be easily 
separated. The angular part of the solution is given in terms of spherical harmonics, and 
the radial part is expressed in terms of modified spherical Bessel functions. For the first 
harmonic, or equivalently, the monopole term in a multipole expansion, and discarding the 
solution growing exponentially, one then finds 

e -m D r 

A = a , Ai = . (4.9) 

r 

which corresponds to a screened radial electric field. 

The non-Abelian Kubo equations have been studied in the literature using different 



Ansatze p3 |. The static case for the SU(2) group with spherical symmetry was considered 
in Ref. [15|. Two particular solutions were found there. The first one corresponds to the 
Yang-Mills vacuum. The second one can actually be constructed by taking the Abelian 
solution in a specific direction of isospin space. Thus 

e -m D r 

A a =f a a — — , A* = 0, (4.10) 
which describes a screened color electric field. Pure non-Abelian solutions were also studied 



in Ref. [jig] , but it was shown there that they tend asymptotically to either the Yang-Mills 



vacuum or to (|4.10| ). 



Let us discuss now the magnetic Kubo equations. In the Abelian case, it is possible to 
find solutions to the Eq. (|2.29|) , just by writing the magnetic field in terms of a "magnetic" 
potential, B = — VA . This is possible since in this case V x B = 0, as we have already 
mentioned. Therefore, in terms of Ao, the field equation becomes 

(V 2 -myi (x) = 0, (4.11) 



which is exactly the same equation as ( |4.8p , and therefore, has the same solutions. In 
particular, with spherical symmetry, the monopole solution is 
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e -m D r 

A = a , A = , (4.12) 

r 

which gives the screened radial magnetic field 

B = -a ^ (I + m D r) e-™^ (4.13) 

Let us consider now the non-Abelian case, for the group £0(3) = S , f/(2)/Z 2 . In this case 
it is possible to construct easily solutions from the Abelian ones. It is enough to consider 
vanishing vector gauge field configurations in all except one specific direction in color space. 
Then the dual non-Abelian Kubo equations reduce to the dual Abelian ones. For example, 
one solution is given by taking in the third direction of isospin space 

B 3 = -«o^ (1 + rh D r) e~^ r , B x = B 2 = , (4.14) 

which describes a screened color magnetic field. 



V. DISCUSSION 

The purpose of this article has been finding a classical mechanism which generates the 
thermal screening of magnetic fields. It has been shown that at finite T, and in the presence 
of magnetic charges, magnetic fields are screened. This effect can be easily understood in the 
Abelian theory. The duality symmetry of electromagnetism allows to derive the magnetic 
screening from the electric one, without further complications. A non-Abelian theory is 
not symmetric under a duality transformation, but still it is possible to show that also 
color magnetic fields can be screened. In both cases, the effects of magnetic screening are 
proportional to the inverse of the gauge coupling constant. 

The question which remains to be answered is whether the polarization effects that have 
been described in this article could be relevant in realistic quantum gauge theories, such 
as QED or QCD. The question can be actually reduced to answer if there are magnetic 
monopoles in those theories. 

It seems obvious that QED does not possess magnetic monopoles. Thus, it is not expected 
that the mechanism described in this article takes place in QED. A pure Abelian gauge theory 
does not suffer from IR problems. The IR divergences of QED that arise in the magnetic 



sector of the theory at finite T can actually be cured in the same way than at zero T [24 



In a pure Yang-Mills theory 't Hooft |25|] showed that in the maximal Abelian gauge that 



theory has magnetic monopoles. Actually this occurs in several different Abelian projections 
of the non-Abelian theory. The idea of describing the QCD ground state as a condensate of 
magnetic monopoles has deserved much attention |26| , p7| , since then one could understand 
the confinement of QCD as a dual Meissner effect. Although these ideas are very appealing, 
there is not yet a complete gauge independent analytical proof of this confinement scenario, 
except for the case of some supersymmetric non-Abelian models However, lattice 

computations have shown the magnetic monopole dominance for the string tension of QCD 
at zero temperature p9 |. 
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A natural expectation is that if magnetic monopoles explain the confinement of QCD 
because they condensate at zero temperature, they should also play an important role in 
explaining the magnetic screening at finite temperature. I have shown that this screening 
would occur by considering a plasma of magnetic charges. Let me stress that the real 
situation corresponding to the real QCD plasma would be much more complicated than 
the simple models that I have described, as then both electric and magnetic polarization 
phenomena should occur simultaneously. There is also a mismatch between the magnetic 
mass of order T/g that has been derived here, and an expected one of order g 2 T. Our result 
has been obtained on the assumption that the equilibrium density of monopoles is of order 
T 3 . However this does not have to be so necessarily. To recover the expected magnetic mass 
of order g 2 T, the equilibrium density of magnetic monopoles should be of order (g 2 T) 3 . This 
would naturally imply that magnetic monopoles are not elementary particles, but that they 
should be dynamically generated. 

Lattice computations could check whether the mechanism of magnetic screening that I 
have described takes place or not in a pure Yang-Mills theory at finite temperature. 
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